
> ################
> #PSRM: Explaining Support for Redistribution: Social Insurance Systems and Fairness
> #
> #Observational Data: Eurobarometer
> #Table C 6
> #
> #Verena Fetscher
> #July 2022
> ####################
> 
> 
> 
> ####################
> # Load data
> ####################
> # EB
> load("EB_benefits_17.Rda")

> ##################
> # Fairness Beliefs and Redistribution Preferences
> ####################
> 
> ##################
> # Table C.6: Linear regression: Support for redistribution, fairness beliefs and benefit concentration
> ####################
> 
> mod1<-lm_robust(redistr~rr+imp_wealthFam+imp_connections+imp_workhard+age+gender+edu_r+leftright+hhsingle+unemp,
+                           data=df_merge[df_merge$incgroup==3,])

> tidy(mod1)
              term     estimate   std.error   statistic       p.value     conf.low
1      (Intercept)  4.429996617 0.185515812  23.8793480 1.485099e-115  4.066245991
2               rr -0.681315266 0.086837508  -7.8458639  5.925306e-15 -0.851582142
3    imp_wealthFam  0.160300543 0.021151348   7.5787387  4.616989e-14  0.118827979
4  imp_connections  0.091697824 0.024314301   3.7713534  1.654889e-04  0.044023492
5     imp_workhard -0.019977401 0.023839725  -0.8379879  4.021041e-01 -0.066721206
6              age  0.004379537 0.001244857   3.5181051  4.410748e-04  0.001938681
7           gender  0.156671302 0.038750261   4.0431032  5.407156e-05  0.080691624
8            edu_r -0.100927273 0.023534324  -4.2885138  1.855319e-05 -0.147072262
9        leftright -0.171688565 0.010308028 -16.6558108  1.188658e-59 -0.191900058
10        hhsingle -0.149098515 0.065490783  -2.2766336  2.287805e-02 -0.277509743
11           unemp  0.265085058 0.657003208   0.4034760  6.866267e-01 -1.023135749
      conf.high   df outcome
1   4.793747243 3009 redistr
2  -0.511048389 3009 redistr
3   0.201773106 3009 redistr
4   0.139372156 3009 redistr
5   0.026766403 3009 redistr
6   0.006820393 3009 redistr
7   0.232650980 3009 redistr
8  -0.054782284 3009 redistr
9  -0.151477071 3009 redistr
10 -0.020687287 3009 redistr
11  1.553305864 3009 redistr

> texreg(list(mod1),stars = c(0.01,0.05,0.1),
+        custom.coef.names = c("Intercept", "Benefit concentration","Wealthy family",
+                              "Connections","Working hard","Age","Female","Education",
+                              "Leftright","Single","Unemployed"))

\begin{table}
\begin{center}
\begin{tabular}{l c}
\hline
 & Model 1 \\
\hline
Intercept             & $4.43^{*}$        \\
                      & $ [ 4.07;  4.79]$ \\
Benefit concentration & $-0.68^{*}$       \\
                      & $ [-0.85; -0.51]$ \\
Wealthy family        & $0.16^{*}$        \\
                      & $ [ 0.12;  0.20]$ \\
Connections           & $0.09^{*}$        \\
                      & $ [ 0.04;  0.14]$ \\
Working hard          & $-0.02$           \\
                      & $ [-0.07;  0.03]$ \\
Age                   & $0.00^{*}$        \\
                      & $ [ 0.00;  0.01]$ \\
Female                & $0.16^{*}$        \\
                      & $ [ 0.08;  0.23]$ \\
Education             & $-0.10^{*}$       \\
                      & $ [-0.15; -0.05]$ \\
Leftright             & $-0.17^{*}$       \\
                      & $ [-0.19; -0.15]$ \\
Single                & $-0.15^{*}$       \\
                      & $ [-0.28; -0.02]$ \\
Unemployed            & $0.27$            \\
                      & $ [-1.02;  1.55]$ \\
\hline
R$^2$                 & $0.17$            \\
Adj. R$^2$            & $0.17$            \\
Num. obs.             & $3020$            \\
RMSE                  & $1.04$            \\
\hline
\multicolumn{2}{l}{\scriptsize{$^*$ 0 outside the confidence interval.}}
\end{tabular}
\caption{Statistical models}
\label{table:coefficients}
\end{center}
\end{table}

> ################
> #PSRM: Explaining Support for Redistribution: Social Insurance Systems and Fairness
> #
> #Observational Data: Eurobarometer
> #Figure C 10
> #
> #Verena Fetscher
> #July 2022
> ####################
> 
> 
> 
> 
> ####################
> # Load data
> ####################
> # EB
> load("EB_benefits_17.Rda")

> ##################
> # Fairness Beliefs and Redistribution Preferences
> ####################
> 
> # Country selection
> df_merge %>%
+   group_by(cntry_short) %>%
+   summarise(wealthFam=mean(imp_wealthFam, na.rm=T))
# A tibble: 13 × 2
   cntry_short wealthFam
   <chr>           <dbl>
 1 AT               3.40
 2 BE               2.71
 3 DE               2.97
 4 DK               2.30
 5 ES               3.01
 6 FI               2.73
 7 FR               2.48
 8 GB               2.81
 9 IE               2.85
10 IT               3.67
11 NL               2.52
12 PT               3.32
13 SE               2.30

> # IMPORTANT FOR GETTING AHEAD IN LIFE - COMING FROM WEALTHY FAMILY (1 Not important at all - 5 Essential)
> 
> ##########################
> #Figure C.10: Benefit concentration and beliefs in the importance of coming from a wealthy family.
> ##########################
> 
> df_merge%>%
+   filter(incgroup==3)%>%
+   ggplot(aes(x = rr, y = imp_wealthFam,label=cntry)) +
+   stat_summary(fun.data = mean_se,aes(color="rich"),alpha=.5)+
+   geom_smooth(method=lm,se=F,color="grey20",na.rm=T,size=.5)+
+   xlab("Benefit concentration")+ ylab("Importance family wealth")+
+   coord_cartesian(ylim=c(2, 4), xlim=c(0,.9))+
+   stat_summary(aes(label=cntry), fun="mean", geom="text", 
+                hjust=0.5, vjust=-.8,size=6)+
+   scale_colour_manual(name="Income group", 
+                       values=c("grey20"),
+                       breaks = c("rich"),
+                       labels = c("rich")
+   )+
+   scale_x_continuous(expand = c(0, 0)) + 
+   scale_y_continuous(expand = c(0, 0)) +
+   theme_bw()+
+   theme(panel.border=element_blank(),axis.line=element_line(),
+         legend.position="none",legend.text = element_text(size = 18),
+         legend.title = element_text(size=18),axis.text=element_text(size=18),
+         axis.title=element_text(size=22)
+   )

> ggsave(file="figure_c10.pdf", height = 5.83, width = 8.27, units = "in")

> ################
> #PSRM: Explaining Support for Redistribution: Social Insurance Systems and Fairness
> #
> #Observational Data: Eurobarometer
> #Figure C 11
> #
> #Verena Fetscher
> #July 2022
> ####################
> 
> 
> 
> ####################
> # Load data
> ####################
> # EB
> load("EB_benefits_17.Rda")

> ##################
> # Fairness Beliefs and Redistribution Preferences
> ####################
> 
> 
> ##########################
> #Figure C.11: Benefit concentration and beliefs in the importance of knowing the right people.
> ##########################
> 
> df_merge%>%
+   filter(incgroup==3)%>%
+   ggplot(aes(x = rr, y = imp_connections,label=cntry)) +
+   stat_summary(fun.data = mean_se,aes(color="rich"),alpha=.5)+
+   geom_smooth(method=lm,se=F,color="grey20",na.rm=T,size=.5)+
+   xlab("Benefit concentration")+ ylab("Importance connections")+
+   coord_cartesian(ylim=c(3, 4.5), xlim=c(0,.9))+
+   stat_summary(aes(label=cntry), fun="mean", geom="text", 
+                hjust=0.5, vjust=-.8,size=6)+
+   scale_colour_manual(name="Income group", 
+                       values=c("grey20"),
+                       breaks = c("rich"),
+                       labels = c("rich")
+   )+
+   scale_x_continuous(expand = c(0, 0)) + 
+   scale_y_continuous(expand = c(0, 0)) +
+   theme_bw()+
+   theme(panel.border=element_blank(),axis.line=element_line(),
+         legend.position="none",legend.text = element_text(size = 18),
+         legend.title = element_text(size=18),axis.text=element_text(size=18),
+         axis.title=element_text(size=22)
+   )

> ggsave(file="figure_c11.pdf", height = 5.83, width = 8.27, units = "in")

> ################
> #PSRM: Explaining Support for Redistribution: Social Insurance Systems and Fairness
> #
> #Observational Data: Eurobarometer
> #Figure C 12
> #
> #Verena Fetscher
> #July 2022
> ####################
> 
> 
> 
> ####################
> # Load data
> ####################
> # EB
> load("EB_benefits_17.Rda")

> ##################
> # Fairness Beliefs and Redistribution Preferences
> ####################
> 
> 
> # IMPORTANT FOR GETTING AHEAD IN LIFE - HARD WORK (1 Not important at all - 5 Essential)
> 
> ##########################
> # Figure C.12: Benefit concentration and beliefs in the importance of working hard.
> ##########################
> 
> 
> df_merge%>%
+   filter(incgroup==3)%>%
+   ggplot(aes(x = rr, y = imp_workhard,label=cntry)) +
+   stat_summary(fun.data = mean_se,aes(color="rich"),alpha=.5)+
+   geom_smooth(method=lm,se=F,color="grey20",na.rm=T,size=.5)+
+   xlab("Benefit concentration")+ ylab("Importance hard work")+
+   coord_cartesian(ylim=c(3.5, 4.5), xlim=c(0,.9))+
+   stat_summary(aes(label=cntry), fun="mean", geom="text", 
+                hjust=0.5, vjust=-.8,size=6)+
+   scale_colour_manual(name="Income group", 
+                       values=c("grey20"),
+                       breaks = c("rich"),
+                       labels = c("rich")
+   )+
+   scale_x_continuous(expand = c(0, 0)) + 
+   scale_y_continuous(expand = c(0, 0)) +
+   theme_bw()+
+   theme(panel.border=element_blank(),axis.line=element_line(),
+         legend.position="none",legend.text = element_text(size = 18),
+         legend.title = element_text(size=18),axis.text=element_text(size=18),
+         axis.title=element_text(size=22)
+   )

> ggsave(file="figure_c12.pdf", height = 5.83, width = 8.27, units = "in")
